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It has been predicted that particles with imaginary mass, called tachyons, would be able to 
travel faster than the speed of light. There has not been any experimental evidence for tachyons 
occurring naturally. Here, we propose how to experimentally simulate Dirac tachyons with trapped 
ions. Quantum measurement on a Dirac particle simulated by a trapped ion causes it to have an 
imaginary mass so that it may travel faster than the effective speed of light. We show that a Dirac 
tachyon must have spinor-motion correlation in order to be superluminal. We also show that it 
exhibits significantly more Klein tunneling than a normal Dirac particle. We provide numerical 
simulations of realistic ion systems and show that our scheme is feasible with current technology. 


I. INTRODUCTION 


In the 1960s, it was theorized that there could be par¬ 
ticles with imaginary mass, called tachyons, which travel 
faster than light m- This is because in order for rela¬ 
tivistic energy E = me?j— {vjc)^ to be real when m 
is imaginary, v must be larger than c. Note that although 
a tachyon is considered to be superluminal, causality may 
still be preserved mu. There has been no experimental 
evidence of tachyons occurring naturally 13 ig. There 
have been a few proposals for engineering tachyon-like 
excitations that travel faster than the effective speed of 
light in the system. These include inverted optical media 
[5] and waveguides [SHII] ■ The motivation for an exper¬ 
imental realization is that it would broaden the range 
of accessible phenomena, and allow the further study of 
physics that would otherwise be unphysical. 

Ion traps have proven to be an ideal setting for exper¬ 
imentally realizing quantum relativistic effects [T^20j . 
Other platforms have been considered as well pTII23j . 
Recent trapped-ion experiments have simulated Dirac 
particles and observed Zitterbewegung and Klein tunnel¬ 
ing [mini- In these experiments, the excellent control 
and read-out capabilities allow one to prepare and mon¬ 
itor the wavepacket dynamics. 


In this paper, we propose a scheme to simulate Dirac 
tachyons |24j with trapped ions and show that their quan¬ 
tum nature distinguishes them from classical tachyons. 
In our scheme, continuous measurement on a Dirac par¬ 
ticle simulated by a trapped ion causes it to have an 
imaginary mass. The Dirac tachyon can then move faster 
than the effective speed of light in the system (see Fig. [I I. 
We perform realistic numerical simulations with example 
experimental numbers and show that this scheme is fea¬ 
sible with current technology. We also describe how to 
measure the relevant observables. 


We also obtain new results regarding the properties of 
Dirac tachyons. We show that spinor-motion correlation 
plays a crucial role in the propagation, i.e., there must 


be spinor-motion correlation in order for a tachyon to be 
superluminal. We also consider the interaction with an 
external potential and find that a tachyon exhibits sig¬ 
nificantly more Klein tunneling than a normal particle. 
Then, we use a spacetime duality to show that a nor¬ 
mal particle scattering off a spatially-varying potential 
is dual to a tachyon scattering off a time-varying poten¬ 
tial. All these features can be simulated with trapped 
ions following our proposal. 

Our paper is structured as follows. In Sec. we in¬ 
troduce the Dirac equation with imaginary mass and de¬ 
scribe our proposed experimental scheme. In Sec. |IIH we 
discuss the velocity of the wavepacket and its relationship 
to entanglement. In Sec. |IV[ we calculate Klein tunneling 
for a Dirac tachyon. In Sec. we provide experimental 
numbers and comment on the scalability of our scheme. 


II. MODEL 

A quantum relativistic particle is described by the 
Dirac equation [25] . In one spatial dimension, it takes 
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FIG. 1. (color online) Propagation of a tachyon wave packet, 
(a) Probability density |'!/)(x,t)p is plotted on color scale. 
Dashed line denotes light cone, (b) Solid line is average po¬ 
sition (x) as a function of t. Dashed line denotes speed of 
light. Initial wavefunction is V’(®) = u+{po) 

with po = 3.5/A and m = 2/(cA). 
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the form idtip = Hip, where [521 

H = cpcjx + mc^az (1) 

with speed of light c, mass m, momentum operator p, 
and Pauli matrices ak- We assume that h = 1. The 
wavefunction ip has both motional and spinor degrees of 
freedom. For the two-component spinor, we work in the 
az basis: |t) and ||). 

This equation was experimentally implemented with 
a trapped ion by letting the ion’s electronic and vibra¬ 
tional degrees of freedom correspond to a Dirac particle’s 
spinor and motional degrees of freedom [151 ITS] , The 
term pa^ can be written in terms of acr± and a^cr±, so 
by exciting both blue and red motional sidebands, one 
obtains Eq. Q. The effective speed of light and mass 
are c = 2ryAf2 and mc^ = fl, where p is the Lamb-Dicke 
parameter, A is the size of the ground-state wavefunc¬ 
tion, n is the strength of the carrier transition, and D is 
the detuning. 

Now suppose It) has a finite lifetime given by a 
linewidth 7 . Then, conditioned on the absence of a decay 
event, the system evolves with a non-Hermitian Hamil¬ 
tonian [57H5D] . 

H = Ho- "^az. (2) 

This is because the null measurement of a decay event 
has back-action on the wavefunction, which is accounted 
for by the non-Hermitian term. (See footnote [3I])- Let¬ 
ting Ho = cpGx, we then have the Dirac equation with 
imaginary mass [24] . 

H = cpax — imc^az, (3) 

where mc^ = 7/4. [Note that the sign of the mass term 
in Eq. ([^ is unimportant since it can be flipped via a uni¬ 
tary transformation.] It is important to condition on the 
absence of decay events; if decay events were included, 
the dynamics would be described by a master equation 
instead of a non-Hermitian Hamiltonian [551 - 151] . Thus, 
quantum measurement allows one to simulate a Dirac 
tachyon with a trapped ion. 

The experimental protocol is as follows. Starting 
from an initial wavepacket, one evolves the ion with 
Ho = cpax, by resonantly exciting blue and red side¬ 
bands, while optically pumping |t) to an auxiliary state. 
After a given amount of time, one turns off Ho and the 
optical pumping, and then measures the population in 
the auxiliary state using the usual fluorescence method. 
If one measures no population in the auxiliary state, then 
there was no decay event and therefore the ion was evolv¬ 
ing solely according to E q. ([^ . This is a probabilistic 
but heralded method [55lI37| : one repeats the protocol 
many times and post-selects on those runs without decay 
events, since those are the ones that simulate Eq. (|^. In 
Sec. [^ we provide experimental numbers and numerical 
simulations with realistic decoherence and show that the 
probability of a successful run is reasonably high. Also, 
it is not necessary to efficiently detect a single photon; 
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FIG. 2. Dispersion relations (energy vs. momentum) for (a) 
normal massive particle and (b) tachyon. Dashed lines denote 
speed of light. 

instead, one only needs to efficiently measure the popu¬ 
lation in the auxiliary state. 

We distinguish our scheme from related work. If the 
second term of Eq. (§ included a complex-conjugation 
operator, as in the Majorana equation [II El [55] , one 
could simulate it by embedding the spinor in a larger 
Hilbert space; however, this does not work for the 
tachyon quantum simulation. Also, a Dirac tachyon is 
different from a Klein-Gordon tachyon HH] because of 
the spinor degree of freedom, which is important to the 
dynamics. In addition, our scheme is different from su- 
perluminality in an absorptive medium with anomalous 
dispersion [SSHU, because the latter is due to pulse re¬ 
shaping instead of imaginary mass and does not have a 
tachyonic dispersion [Fig. [^b)]. 

We now compare the behavior of normal massive par¬ 
ticles [described by Eq. Q] with tachyons [described by 
Eq. (§]. 

III. VELOCITY 

A normal massive particle travels slower than c, be¬ 
cause the group velocity Vg = dE/dp is smaller in magni¬ 
tude than c, where E = -I- m?c'^ is the dispersion 

of Eq. Q, plotted in Eig. [^a). Another way to see this 
is from the equation of motion for the position operator, 
dt{x) = c{ax), which is bounded above by c. 

The special property of a tachyon is that its group 
velocity is larger than c. This is because the dispersion 
of Eq. 0 is E = E\/p^(? — m?c'^, as shown in Fig. [^b). 
Figure]^ is a spacetime plot of a tachyon, showing that 
the wave packet is indeed superluminal. 

Consider the equation of motion of {x) for a tachyon. 
For a non-Hermitian Hamiltonian like Eq. the equa¬ 
tion of motion for an operator includes extra terms mi, 
leading to 

dt{x) = c{ax) - 2m(?{{xaz) - {x){az)). (4) 

Since the first term is bounded by c, in order to have 
|9t(a:)| > c, the quantity in parentheses (correlation of 
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FIG. 3. (color online) (a) Correlation of x and CTz for Fig. Bb). 
(b) Population in |t) (solid line) and |^) (dashed line) for the 
state in Eq. with Pa = 3.5/A and m = 2/(cA). 


X and cr^) must be nonzero. In other words, in order 
for a Dirac tachyon to be superluminal, the spinor and 
motional degrees of freedom must be correlated. Without 
correlations, a Dirac tachyon is actually subluminal. (For 
the case of a pure state, superluminality requires spinor- 
motion entanglement, but in an experiment, the state 
will be mixed.) 

This explains the behavior in Fig. m- The initial 
wavefunction is tp{x) = u+{po), where the 

spinor u+{po) is the positive-energy eigenstate of Eq. ^ 
for p = Pa. Since this is a separable state, the tachyon 
is initially subluminal. Over time, the wavefunction de¬ 
velops correlations [Fig. [^a)] so that it eventually be¬ 
comes superluminal. This oscillatory motion is an ex¬ 
ample of Zitterbewegung [55] and occurs because the 
initial wavefunction contains both positive and negative- 
energy components that interfere with each other. In 
momentum space, ipip) is a Gaussian centered at po, but 
M+(Po) is the positive eigenstate only for p = p^- 

The following entangled state, 

ipix) = J dpe~^ (p-Po) e*P“u_|_(p), (5) 

contains only positive-energy components because the 
spinor u+{p) is the positive eigenstate for each p. Thus, 
this wavepacket travels superluminally with group veloc¬ 
ity Vg{po) and without Zitterbewegung. For motion to the 
right, X and CTz are negatively correlated, i.e., the peak 
of |V’i(a;)p is ahead of \ijj'f{x)\'^, as shown in Fig. rab). 
In the limit of Po ^ me, (xa^) — {x){az) = —mc^p^- 
In contrast, for a normal particle in the state given by 
Eq. m, (xCTz) - {x){az) = 0. 


IV. SCATTERING OFF A POTENTIAL 

Here, we consider the scattering off a repulsive elec¬ 
tric potential that is linear in space, e4>{x) = gx, 
where e is the charge and g > 0. It is well known 
that a normal Dirac particle can tunnel through such 
a barrier and propagate undamped due to the negative- 
energy branch of the dispersion [Fig. [^a)]. An incoming 


wavepacket splits into a reflected component and a tun¬ 
neling component [Fig. Ba)], and the tunneling proba¬ 
bility is Ptunnei = exp(—TTm^c^/p) [ 32 ] • This is known as 
Klein tunneling |43j and is surprising because a nonrel- 
ativistic Schrodinger particle would not tunnel without 
attenuation. In the trapped-ion implementation, a lin¬ 
ear potential can be engineered by using a second ion, as 
demonstrated recently [HHI]. 

Now consider a tachyon scattering off a linear electric 
potential [Fig.^b)], 

idtip = {epax — im(?az + gx)'ijj. ( 6 ) 

Suppose an initial wavepacket moves to the right with 
p = Pa > 0 and E > 0. We work in momentum space 
and write x = idp and ' 0 (p, t) = ^{po = p + gt, t) [T5] : 

= [c(Po - 9t)<Jx - ime^az]^. (7) 

Thus, the motion of the wavepacket is equivalent to a 
Landau-Zener process with a magnetic held c{po — gt) 
that decreases linearly in time [44]. The energy levels 
of this process are given by the dispersion in Fig. [^b). 
At t = 0, the system starts at p = Pa on the positive- 
energy branch, p is then ramped to negative values, and 
the question is how the hnal population is distributed 
between the branches. The population in the positive- 
energy branch corresponds to the (left-moving) rehected 
component, while the population in the negative-energy 
branch corresponds to the (right-moving) tunneling com¬ 
ponent. [Note that the eigenvalues for p € {—me, me) are 
complex and not plotted in Fig. Bb)]. 

However, this is a non-Hermitian Landau-Zener pro¬ 
cess due to the imaginary held in Eq. Q. By extending 
Zener’s original solution [ 33 ] to the non-Hermitian case, 
we hnd that the tunneling probability is 



This result differs from that for a normal particle: the 
tunneling probability for a tachyon is larger and is al¬ 
ways at least 1/2. For example. Fig. shows that a lot 
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X (units of A) 


FIG. 4. (color online) Probability density \'tp{x,t)\^ scat¬ 
tering off a linear electric potential for (a) normal massive 
particle and (b) tachyon. For both plots, m = l/(cA) and 
g = 2c/A^ 
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more population tunnels through for a tachyon than a 
normal particle. Physically, this is because the tachyonic 
dispersion does not have an energy gap, while the normal 
dispersion does (Fig. [^. 

There is a useful duality between normal Dirac parti¬ 
cles and tachyons. Suppose is a solution to the 

Dirac equation for a normal particle with electric po¬ 
tential (^(cc) and vector potential ^(a;) [121 US]. One 
can show that tj}{x',t') = U~'^tp{x = t',t = x') is a 
solution to the Dirac equation for a tachyon with elec¬ 
tric potential (/''(t') = —A{x = t') and vector potential 
A'{t') = —(j){x = t'), where U = is a unitary 

transformation. Thus, a normal particle scattering off a 
spatially-varying electric potential is dual to a tachyon 
scattering off a time-varying vector potential. 


(a) 




FIG. 5. (color online) (a) Level scheme for (b) 

Nnmerical simnlation of trapped-ion implementation for dif¬ 
ferent rates of optical-pnmping errors. 


V. EXPERIMENTAL CONSIDERATIONS 

A. Probabilities 

The tachyon protocol is probabilistic (i.e., conditioned 
on no decay event), so it is important to estimate the 
probability of a successful experimental run. To do this, 
we note that the average number of decay events during 
time t is 

+ -ft 

^ = -2-"T’ 

where we have used the fact that the eigenstates of 
Eq. § satisfy (a^) = 0. The number of decay events 
is a Poissonian random variable. Thus, the probability 
of a successful run (zero decay events during time t) is: 

^success — e ^ ~ e ^ . (fl^) 

(To calculate ^success exactly, one would use the instan¬ 
taneous value of (iTz), but we have found that setting 
(cTz) = 0 in Eq. results in a very good estimate.) 

We recall from Eq. ® that 7 = 4toc^. We also note 
that m is in units of l/fcA), while t is in units of A/c 
(since h = 1): 



where m',t' are dimensionless. Thus, 

^«2mY, Psuccess«e-2™'‘'. (12) 

As an example, consider Fig. [T| where m! = 2 . The 
particle is superluminal by t' = 0.5, which corresponds 
to Psnccess = 6 “^ « 0.14. A time of t' = 1 corresponds to 
^success = « 0.018. These probabilities are reason¬ 

ably high enough for the scheme to be feasible. Psuccess 
can be easily increased by decreasing m'. Also, we note 
from Eq. Q that if the initial wavefunction is suitably 
entangled ]45j . the initial velocity can be larger than c, 
so that superluminality occurs earlier. For the case of an 
external potential, one can observe Klein tunneling with 
^success « 0.03 by using a large g. 


B. Example numbers 

We provide example numbers for See the 

level scheme in Fig. [^a). To implement the first term 
of Eq. one drives the red and blue sidebands of the 
|4,)-|t) transition (either directly or via a P 1/2 state). 
To implement the second term of Eq. (|^, one optically 
pumps It) to the auxiliary state |a) via the IP 3 / 2 , P = 2 ) 
state. Due to dipole-selection rules, |P 3 / 2 ,P = 2) de¬ 
cays into I a) instead of |t). Also, the branching ratio 
of 1 ^ 3 / 2 ) E = 2 ) to It) is less than 0.002 US] ST], so these 
dephasing errors can be neglected [Fig.j^b)] . To measure 
the population in |a), one excites |a) to^Pi/ 2 , P = 0) and 
looks for fluorescence. To measure the ion position over 
time, one uses the procedure described in Refs. [HI El- 

Setting T] = 0.05, A = 3.4 nm, and D = 27r x 100 kHz 
means an effective speed of light c = 2 x 10 “"^ m/s and 
time scale A/c = 16 /xs. For optical pumping 7 = 27r x 
80 kHz, the mass is m = 2/(cA). These numbers are 
readily accessible in current experiments. 


C. Numerical simulatious 

We have performed numerical simulations of the 
tachyon protocol including realistic decoherence present 
in ion traps. The carrier term is eliminated by using 
two counterpropagating laser fields for red and blue side¬ 
bands, Ho = — a+e~'^^^)[sm (j)+r] cos + 

a'^e®'^*)], where v is the trap frequency and the laser 
phase (p is chosen such that sin^ = 0. For the tachyon 
protocol, the intrinsic source of decoherence is optical 
pumping errors, i.e., when |P 3 / 2 ,P = 2 ) decays back to 
It), which leads to dephasing with rate 7 ^. Figure [^b) 
shows that the optical pumping errors are negligible for 
171 Yb+, since 7 ^ = O.OO 27 . We have also found that 
since 7 ^ is much larger than other sources of decoher¬ 
ence (motional heating, spontaneous emission, and laser 
dephasing) in current experiments, the latter have even 
less effect. Thus, the tachyonic behaviour is preserved 
with realistic decoherence. 
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D. Dependence on detector efficiency 

Since our scheme is conditioned on the absence of a 
decay event, one must reliably measure whether a decay 
event occurred. In a typical ion trap experiment, the 
efficiency of detecting a single photon is very low. Fortu¬ 
nately, the scheme described in Sec. |VB| does not rely on 
detecting a single photon. The key is to have |t) decay 
into I a) instead of ||). To determine whether a decay 
event occurred, one measures the population in |a) by 
continuously exciting it to 1^1/2, F = 0). If the atom is 
in I a), it will scatter many photons; otherwise, it will not 
scatter any photons. This way, one can efficiently read 
out the population in |a) and thereby determine whether 
a decay event occurred, even though the single-photon 
efficiency is low. In other words, instead of detecting the 
single photon created during the decay from |t), one can 
just measure the many photons that are scattered when 
the atom is in |a). 

How high should the readout fidelity be to observe 
tachyonic behavior? If there is a readout error, it is sim¬ 
ilar to an optical pumping error discussed in Sec. |V C[ 
i.e., one measures that there was not a decay event even 
though there really was. According to Fig. a readout 
fidelity of 99.5% is sufficient for observing that the group 
velocity exceeds the effective c. For comparison, a read¬ 
out fidelity of 99.99% has been achieved with a trapped 
ion [48] . 


E. Measurements 

Spinor-motion correlation is important to tachyonic 
dynamics [Eq. so we describe a protocol to directly 
measure it. It is well known that single-qubit observ¬ 
ables {(Jx,y,z) can be measured by resonance fluorescence 
of a cycling transition. This procedure can also be used 
to measure vibrational degrees of freedom by mapping 
them onto spinor ones. The expectation value {xa^} is 
measured by observing (cr^) as follows. A combination of 
red and blue sidebands, Hr = riQ{aa+e^^^ + 
and Hb = -I-a(T_e“®‘^*’), implements a state- 

dependent displacement when the phases are properly 
tuned. The corresponding dynamics, U = acts 

on az according to Ua^W = cos{kx)(Tz + sm{kx)ay = O. 
For k{x) <C 1, this can be approximated by Icr^ -|- kxay. 
Therefore, the derivative of the expectation value of O, 
dk{0), corresponds to {xoy), which can be transformed 
into {xUz) with a rotation of the spinor. 


F. Scalability 

In this paper, we discussed how to simulate a sin¬ 
gle Dirac tachyon using a single trapped ion. Here, we 
consider the possibility of generalizing the scheme to N 


ions. Suppose there is a one-dimensional array of N 
trapping potentials (as in a segmented trap), each with 
an ion inside. Due to the normal-mode structure from 
the Coulomb interaction, one applies lasers to induce a 
mode-mode coupling. Further lasers drive blue and red 
sidebands, and one conditions on the absence of a de¬ 
cay event. In this way, one implements a non-Hermitian 
many-body system. There are two questions. First, how 
large could A^ be in a realistic experiment? Second, how 
large does N have to be in order to be difficult to simulate 
on a classical computer? 

For N ions, Psuccess ~ ‘ ^. Suppose we set m' = 

0.25 and t' = 1. Then for N = 5, Psuccess = 0.08. This is 
still reasonably high to be experimentally realistic. 

Now, to simulate such a system on a classical computer 
would require representing the wavefunction on a Hilbert 
space. The spin-1/2 degree of freedom has dimension 2. 
For the vibrational degree of freedom, the Hilbert space 
has to be truncated at some number of phonons, but 
there should be enough phonons to accurately represent 
the wavefunction. For example, a wave packet with a 
moderate value of p would require 30 phonons. Thus, 
the Hilbert space for each ion has dimension 60, and the 
total Hilbert space has dimension 60^. For N = 5, the 
dimension is 8 x 10®, which is difficult to simulate on 
a classical computer. (This is equivalent to having 30 
spin-1/2 particles.) Thus, it is possible to reach a regime 
where N is large enough to be difficult to simulate clas¬ 
sically. This is due to the large number of phonons that 
must be tracked. 


VI. CONCLUSION 

We have proposed a quantum simulation of tachyons 
with current ion-trap technology. This ability to ex¬ 
perimentally observe tachyonic physics opens the door 
to many possibilities. For example, one can study how 
imaginary mass affects bound states in a confining po¬ 
tential [IH] or the Dirac oscillator EO). One can also 
consider a tachyon interacting with another tachyon or 
a normal particle |5]. By including more ions, one can 
study many-body physics iMiisziin]- Finally, tachyon 
physics in conjunction with metamaterials may lead to 
new applications in manipulating the propagation of light 

m- 
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